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Abstract 

We numerically study the magnetization and the dispersion relation 
of a frustrated quantum spin system. Our method, which is named the 
stochastic state selection method, is a kind of Monte Carlo method to give 
eigenstates of the system through statistical averaging processes. 

Using the stochastic state selection method with some constraints, we 
make a successful study of the spin-1/2 Heisenberg antiferromagnet on 
a 48-site triangular lattice. We calculate the sublattice magnetization 
and the static structure function in the ground state. Our result on the 
sublattice magnetization is consistent with the value given by the linear 
spin wave theory. This adds an evidence for the analysis based on the 
spontaneous symmetry breaking of the semi-classical Neel order in the 
ground state. 

We also evaluate the low-lying one magnon spectra of the model with 
all wave vectors available on a 48-site triangular lattice. We find that 
at the ordering wave vector there is a Goldstone mode, which is in good 
agreement with the result from the spin wave analysis. The magnon 
spectra with other wave vectors, however, are quite different from results 
obtained by the linear spin wave theory. We observe a flat dispersion rela- 
tion with a strong downward renormalization. Our results are compatible 
with those recently reported in the series expansion study and in the order 
1/S calculation of the spin wave analysis. 

1 Introduction 

Since the pioneering work by Anderson and FazekaspQ |2J, the spin-1/2 Heisen- 
berg antiferromagnet on a triangular lattice has been extensively investigated 
as a promising candidate to realize a spin-liquid ground state induced by geo- 
metric frustration and quantum fluctuations. Yet, in spite of a large amount of 
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theoretical and experimental works, we do not have any unified picture for this 
system. 

On the theoretical side, most of the numerical studies carried out over the 
past decade with a variety of different techniques do not support that the sug- 
gested spin-liquid ground state is realized in this model. Instead they pro- 
vide evidences to indicate the ground state with the three-sublattice order 
where the average direction of neighboring spins differs by a 120° angle[3l [4] 
El El [3 El El Uni HH HH H31 HH US]- Then the linear spin wave theory 
(LSWT)[inj Q21 QH US [13 well describes numerical results calculated on lat- 
tices with finite sizes. 

On the experimental side, several novel materials with triangular structures 
have been investigated recently. One of these materials is CS2CUCI4 [21} [22], 
which is supposed to reduce to the one-dimensional spin- 1/2 quantum Heisen- 
berg antiferromagnet because of its anisotropy 23, 24. Other interesting materi- 
als are /c-(BEDT-TTF) 2 Cu 2 (CN) 3 [2S] and EtMe 3 Sb[Pb(dmit) 2 ]2 HUGH], which 
are considered to be close to the Heisenberg antiferromagnet on an isotropic tri- 
angular lattice. These materials, however, do not show any magnetic long-range 
order down to the quite low temperature compared with the exchange interac- 
tions. 

Through further studies motivated by these experiments, theorists have 
found that fundamental properties on a triangular lattice are quite different 
from those on a square lattice, while antiferromagnets on both lattices have the 
semi-classical long-range orders. The dispersion relation is one of the properties 
that have been investigated to compare systems with different geometries. Re- 
cently the series expansion study [28, 29[ and the full 1/S calculation of the spin 
wave theory(30l EU [32] on this relation show that on a triangular lattice one 
sees a downward renormalization of the higher energy spectra, while on a square 
lattice one sees an upward renormalization. The former authors also point out 
that the roton minimum is present in the relatively flat region of the dispersion 
relation on the triangular lattice. These features are quite different from the 
predictions of the LSWT. 

In these somewhat confusing situations one needs unbiased numerical studies 
which do not depend on any special physical assumption. The stochastic state 
selection (SSS) method, a new type of Monte Carlo method which we have 
developed in several years [33] [3H ESI [36] [37] (38] [39] , has such a good propertj0- 
One can therefore employ the method to evaluate any quantity in any system 
one wants to investigate. 

In the algorithm of the SSS method we consider the full Hilbert space of the 
system and stochastically reduce it to relatively small one so that we can manage 
it in numerical studies. More concretely, we use a set of stochastic variables 
which are as many as basis states of the whole vector space under consideration, 
but most of these variables are valued to be zero. Then any inner product related 
to an arbitrary operator is calculable using the survived non-zero stochastic 

1 AU one assumes in numerical calculations is that the basis states have several general 
symmetries attributed to geometrical symmetries of the lattice on which the system lives. 
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variables. Statistical averaging processes guarantee in a mathematically justified 
manner that the result becomes the correct value of the inner product. It is 
found that several constraints on the set of stochastic variables are helpful to 
obtain better results with less statistical errors. Using this constrained SSS 
method we started our numerical study on the spin-1/2 quantum Heisenberg 
antiferromagnet on a 48-site triangular lattice. We have estimated lowest energy 
eigenvalues of the model for each sectors with < S z < 4, where S z denotes the 
z component of the total spin S of the system [55]. 

In this paper we make a further investigation of the model by means of the 
constrained SSS method with two new applications. One of them is to accu- 
rately calculate expectation values of operators which contain many off-diagonal 
elements in their representations. By evaluating the sublattice magnetization 
and the static structure function we demonstrate that it is possible to obtain 
accurate knowledge of the ground state in this method. It should be noted that 
in the usual quantum Monte Carlo method these physical quantites are not easy 
to calculate even for non-frustrated systems. Another is an extension to employ 
a new set of basis states with which complex coefficients are inevitable in an 
expansion of an arbitrary state. Using this set of basis states in the constrained 
SSS method we successfully calculate low-lying one magnon spectra with non- 
zero wave vectors. It should also be noted that even for non- frustrated systems 
such as the quantum Heisenberg antiferromagnet on a square lattice we cannot 
do without complex numbers in calculations with non-zero wave vectors. Our 
study in this paper performed by means of the constrained SSS method gives 
reliable results calculated from the first principle. We see that our results are 
compatible with those in refs.[2H EH 1301 El [32]. It therefore supports the real- 
ization of an ordered ground state in the model. At the same time, however, it 
adds an evidence that dynamical properties of the system are not described by 
the LSWT. 

The plan of this paper is as follows. In section [2] we make brief descriptions 
of the model and the method. Subscction l2.1l is to define the Hamiltonian of the 
model we study. In addition we comment on the power method. An operator 
related to the Hamiltonian is introduced here so that we can obtain the lowest 
eigenvalue of the Hamiltonian using the power method. In subsection 12.21 the 
SSS method is shortly reviewed. Here an introduction of a two- valued stochastic 
variable that might be zero is essential. Then the constrained SSS method is 
explained. In section [3] we present our numerical results. Subsection 13.11 is to 
show our definition and results on the sublattice magnetization and the static 
structure function. These results are compared with those reported on smaller 
lattices. On the sublattice magnetization we also make a comparison with results 
calculated by the LSWT. In subsection 13.21 we show our results on the one 
magnon spectra. For nine independent states in the Brillouin zone of a 48-site 
triangular lattice we present estimated values of energy eigenvalues. Then the 
one-magnon dispersion relation is studied in comparison with results by the 
LSWT. Section [4] is devoted to summary. Here we also count up characteristic 
features of our method. Finally an appendix is given to explain our sets of basis 
states in detail. 
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2 Model and method 



2.1 Model and power method 

The Hamiltonian of the model on the TV-site lattice is 

i? = J^I,-l,, (1) 

(ij) 

where Si denotes the spin- 1/2 operator on the z-th site and the sum runs over 
all Nb(= 3N) bonds of the iV-site lattice. The coupling J is set to 1 throughout 
this paper. In our calculations we employ the power method. The SSS method 
is used to calculate expectation values of powers of an operator Q, 

Q = ll- H, (2) 

where / denotes the identity operator and I is a positive number. For each 
quantum number we choose one value of / which ensures that the lowest energy 
eigenvalue of the system corresponds to the eigenvalue of Q whose absolute value 
is the largest. Detailed explanations are given in a previous paper [39]. 
In the power method we can estimate the energy eigenvalue E by 

Q E = I - E = lim (ip A | Q m | Va>/W>a I Q m ' X | Va> , (3) 
where | ^a) denotes a state which approximates the exact eigenstate | ^e) of 
I V>a) =| ^e)c+ I Qs (N«|c|), (4) 

since 

lim QT I Va) oc I ^) (5) 

m^oo 

holds. 

We estimate an expectation value of an arbitrary operator O in the eigenstate 
| Ve) using 

(6) = lim (V> A I OQ m | ^ a )/(Va I QT I V>a) • (6) 

m— too 

Note that in this calculation a contamination appears at s/c order because 

(d) = (ite\6\fa) + (ii> E \d\0s/c. (7) 

2.2 Constrained SSS method 

The stochastic state selection is realized by a number of stochastic variables. 
Let us expand a normalized state | tp) in an TVy-dimensional vector space by a 
basis {| j)}, 

N v 
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Then for each j we generate a stochastic variable r/j following to the on-off 
probability function, 

P .( v ) = ±5(r) - a,) + (1 - -)<%), - = min (l, ^f] , (9) 

with a positive parameter e which is common to all Pj{if) {j — 1, 2, • • • , Ny). 
Note that © implies r/j = cij (> 1) or r]j = 0, statistical averages being ({r]j)) = 1 
and {{r]j}} = ctj. In numerical studies we can neglect basis states whose i]j=0. 
The parameter e therefore controls the size of the partial vector space we should 
consider in each sampling. Then calculations in the full vector space are to be 
recovered through statistical averaging processes. 
A random choice operator M is then defined by 

M = Y J \j)m(j\- (io) 

i=i 

Using this M we obtain a state | ip), 

N v 

\^)^M\4,)=Y,\3) C W, (11) 
i=i 

which has less non-zero elements than | -0) has. It should be kept in mind that 
with any operator O the statistical average (( | O | ^2))) is exactly equal to 
the expectation value (tpi \ O | ^2) because ((r]j)) = 1 for all j. 

In calculations using the power method we generate independent random 
choice operator^] M^ 2 \ ■••to obtain an improved state | ip^) f rom 

I Va)- Namely, we define | V'a ) =1 V'a) and 

I } ) = QM {n ~ 1] ■ ■ ■ QM {1) I Va) (n>2). (12) 

In the constrained SSS method, we introduce K (K = 1,2, • • •) constraints 
which can be fulfilled by making K r/j s depend on other N — K independently 
generated stochastic variables for each random choice operator. A detailed 
description on this improvement is given in our previous paper [59"] together 
with numerical studies with a K = 1 case. In the same manner as ref.[3H], we 
use the state | (p& ) instead of the above mentioned | ipj? )• Here | (f>^) is 
defined by | (f)^) =\ tpx) an d 

\^) = QM^...QM^ |Va) (n>2), (13) 

with random choice operator^] M^, , • • •, each of which is generated under 
a constraint 

(V>A I Mj k) I <t>f ) = (V» A I <t>f) (fc-1,2,---). (14) 

2 Superscripts for random choice operators are added here to emphasize that for all j the 
stochastic variable rjj in AfCO is generated independently to the one in JwO with k ^ I. 

3 Here we add a subscript in order to notify that each random choice operator is now 
determined to satisfy a given constraint. 
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We can evaluate an expectation value of an operator O by a statistical av- 
erage of the ratio 

(m) _ (ipA I 6 I jg) 
VpA. I </>A > 

with large enough value of m. Actually, in order to save the computing time, 
either of the following two ratios is employed in our numerical study instead of 
(USD , 

?(m) = (ifo I QMi m) | ^ 

(v-a i a 

or 



R{m) ^ (v» A I m' c o | ^x_0 (17) 

(^A|0i m) ) ' 

where M' c denotes another random choice operator, which is generated indepen- 
dently to M c (1) with the constraint (ip A \ M' c \ ip^) = {ip A \ ip A 1] ). 

It should be noted that the method is applicable for a state which has com- 
plex coefficients as well as a state with real coefficients. For a complex coefficient 
Cj = \cj\e l6j we obtain, after the stochastic selection, CjT]j — or Cj-qj = ee j if 
\cj | < e. With complex coefficients we newly have two technical difficulties. One 
of them is that we need double computer memory to treat complex numbers. 
Another is that we cannot use the symmetries of rotation and reflection which 
reduce the number of independent basis states. 



3 Results 

3.1 Sublattice magnetization and static structure function 

In this subsection we present our results of the sublattice magnetization and 
of the static structure function in the ground state, which we evaluate on the 
48-site lattice using the ratio defined by (fTT)) . We start from an approximate 
state constructed in the sector with S z — since we want evaluate them in the 
ground state of the system. Detailed explanations of the way to calculate this 
approximate state are given in a previous paper [39], 

Since the squared sublattice magnetization, which we denote by Mjy for the 
system with N sites, is defined in ref.jlD] as 



^2 ^2 ^2 

S a + St> + Sr 



6 U ' L ' 

the operator O in ([15)) for should be 



^2 ->2 ->2 

3 ( S A + 5 R + S r 



° = N(N + 6) > (19) 



G 



where A, B and C are sublattices of the TV-site triangular lattice and 



l x = & ( X = A ' B > C ) ■ ( 2 °) 

iex 

Our estimation of the sublattice magnetization with TV = 48 is = 0.380 ± 
0.05. Fig. [T]plots this result together with those reported on TV = 21, 27 and 36 
lattices [3D). As is shown by calculations up to the 36-site lattice, the finite-size 
analysis on the magnetization suggests a finite value of the magnetization in 
the thermodynamic limit. We see that our result is consistent with those for 
smaller lattices, which supports the scaling predicted by the LSWT, 



where ao = 0.2387, bo — 1.215 and Co = — 1.5[41). The finiteness of the sublattice 
magnetization, > 0, is also observed in Fig. [T] 

Next we show results of the static structure function S(k) in the ground 
state, for which the operator O is defined by 



6 = S(-k) ■ S(k) (22) 

with 

S(k) = - 

'N 



Possible values of |fc|/fco on the 48-site lattice, where fco=47r/3, are obtained 
from the Brillouin zone shown in Fig. H They are 1/4, y/3/4, 1/2, ^7/4, 3/4, 
v3/2, VT3/4 and 1. Numerical results with a normalization factor 8/(iV + 6) 
are plotted in Fig. [3] as a function of |fc|/fco- In order to make a comparison 
between our data and previously reported ones for a smaller TV, we also plot the 
structure function calculated on an N = 36 lattice[5] in Fig.[3J We observe that 
in the full range of |fc|/&o our TV = 48 results lie slightly below those calculated 
on the 36-site lattice. 



3.2 One magnon spectra 

In this subsection we discuss the one magnon spectra with S z = 1- Detailed 
descriptions of wave vectors k and a set of basis states {| j(k))} we use here are 
given in the appendix. 

For each value of k = k x e x + k y e y on the 48-site lattice, we start from an 
approximate state | V^a^)) expanded by a set of basis states {| j(k)}} with 
complex coefficients, 

jVv(fc) 

l*A®>= E \J(k))c m , (24) 
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which is used as | -0a) in fT3"j) . (ITU) and so on. The state | ^\(k)) is prepared 
following descriptions in a previous paper [55]. In | 'J'a^)) we constructed for 
each k, the number of basis states with non-zero Cj/j:y which we denote by 
ATbasis) ranges from 270 million to 1.6 billion. We evaluate an upper bound of 
the lowest energy eigenvalue E{k x ,k v ), using the ratio Rq^j in (fH))) with the 
operator O = Q = 11 — H. The results are presented in Table [TJ together with 
values of / and A^asis- 



point 


(k x , ky) 


Z 




E(k x , ky) 





(0,0) 


5.2 


1.6 x 10 y 


-26.327 ±0.023 


Qi 


(vr/3,0) 


5.65 


4.0 x 10 s 


-25.413 ±0.048 


Q 2 


(27r/3,0) 


5.8 


3.7 x 10 s 


-25.202 ±0.065 


Q 3 


(tt.O) 


5.8 


3.5 x 10 s 


-25.241 ±0.061 


Q 


(4tt/3,0) 


5.2 


9.6 x 10 8 


-26.131 ±0.030 


Bi 


(7r/2,ir/(2V3)) 


5.8 


3.2 x 10 8 


-25.248 ±0.076 


B 


(7r,7r/V3) 


5.8 


2.7 x 10 8 


-25.419 ±0.062 


Di 


(57r/6,7r/(2V3)) 


5.6 


3.6 x 10 8 


-25.288 ±0.040 


D 


(77r/6,7r/(2V3)) 


5.65 


3.7 x 10 8 


-25.416 ±0.050 



Table 1: Results for one magnon spectra for values of k at the points in the 
Brillouin zone shown in Fig. [5] Here iVbasis is the number of basis states we 
actually employ to construct the initial approximate state | \&A(fe))- The value 
of Z in Q defined by ([2]) is determined so that it gives the lowest eigenvalue of 
ff[3"9"] with the given k. A typical value of e, the parameter in ©, is 0.002, 
with which ~ 1.3% of the basis states survive after each stochastic selection. 
Depending on k, values of m we employ to calculate the energy upper bound 
lie between 15 and 20, and numbers of samples to calculate the average and the 
statistical error ranges between 100 and 200. 

In Fig. 2] we plot a gap energy of one magnon with a wave vector (k x ,k y ), 
which is calculated from data in Table [JJ by 

A(k x ,k y ) = ~E(k x ,k y ) - £7(0, 0) , (25) 

along the path OQDBO in Fig. [2] We also plot the dispersion relation calculated 
by the LSWT(2"0"] in Fig. H) We see that our data is consistent with the gapless 
state at the ordering wave vector, which is expected by the full breaking of the 
SU(2) symmetry. The gap energy at the central region of the Brillouin zone 
is much smaller than that calculated by the LSWT analysis. This downward 
rcnormalization is in contrast with the upward renormalization found in the 
square lattice. In the same central region we also find that the spectra is flat or 
dispersionless. 
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4 Summary 



In this paper we have reported further developments on numerical calculations 
by means of the stochastic state selection (SSS) method, which has the following 
good properties. 

• With this method we can study large systems compared to the exact 
diagonalizationJB] . 

• This method presents a solution for the so-called sign problem. Note that 
the ordinary quantum Monte Carlo method, which has been quite powerful 
to study spin systems on bipartite lattices [HI H3J H3J, does not work in 
most frustrated systems. 

• One does not need any physical assumption. This is in contrast to the sit- 
uation that recent variational approaches are based on some assumptions 
which are specific to the system to be studied. For example, the coupled- 
cluster methodJTJJ Q3] is based on an eigenstate with a special form. In 
the density matrix renormalization group method [TBI EH SSI SSI Wf\ . 
on the other hand, one assumes that the ground state is represented by 
the matrix product state. Another example is the stochastic reconfigura- 
tion method[lS] whose origin is assigned to the fixed-node Monte Carlo 
method 49 , 50 . In these methods an effective Hamiltonian is used instead 
of the exact one. 

• In the SSS method it is easy to treat complex coefficients in the expansion 
of states. This is impossible in the ordinary Monte Carlo method. 

Using this method we numerically study the spin-1/2 quantum Heisenberg 
antiferromagnet on the 48-site triangular lattice. We evaluate the sublattice 
magnetization and the static structure function in the ground state of the sys- 
tem. We also calculate the one magnon spectra with non-zero wave vectors. 

Our results on the magnetization and the structure function show that the 
SSS method is quite effective to estimate these quantities. The accuracy of the 
data, which depends on how many samples are available to calculate statistical 
averages, is satisfying within our computing resources. We observe that the re- 
sults we obtain are consistent with the reported ones on smaller lattices. There- 
fore we see that the linear spin wave theory describes the results well so that 
an analysis based on the spontaneous symmetry breaking on the semi-classical 
Neel order ground state is acceptable. 

On the magnon spectra, we make sure that the one magnon energy at the 
ordering wave vector is consistent with the Goldstone mode, as is expected in the 
spin wave analysis. For other values of the wave vectors, however, we observe 
flat spectra with strong downward renormalization, which are quite different 
from the spectra predicted by the LSWT. We find our results are qualitatively 
consistent with results by the series expansion study[28j [29] and by the order 
1/5 calculation of the spin wave analysis [501 15T1 152"] . 
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Finally let us again emphasize that we can get physical insight for the frus- 
trated quantum spin systems from numerical results by the SSS method because 
the method is well developed to calculate various observables such as energy, 
the magnetization and the magnon spectra. Especially it is notable that using 
this method we can evaluate the magnon spectra. While neither the ordinary 
Monte Carlo method nor the variational Monte Carlo method is applicable to 
calculate this quantity, we can use the SSS method in the same manner as the 
one used to evaluate the ground state energy. 



Appendix 

In this appendix we summarize our definition of wave vectors k and two sets 
of basis states {| j)} and {| j(k))} employed in this paper. The former set {| j)}, 
which is introduced for the 48-site lattice in the previous work with k = 39 , is 
used to calculate the sublattice magnetization and the static structure function. 
While the latter set {| j(k))} is newly introduced here in calculations of the one 
magnon spectra with non-zero wave vectors. 

Let us denote two basis vectors of the triangular lattice by u\ and ui , which 
are described by unit vectors e x and e y in x- and y-directions, 

ui = e x , u 2 = -e x + —e y . (26) 

The wave vectors are defined by 

k = + -^§ 2 = k x e x + k y e y , (K a = 0, 1, 2, • • ■ (a = 1, 2)) , (27) 

where g a (a = 1,2 ) are two basis vectors of the reciprocal lattice for which 

u a ■ g b = 2Trd ab (a = 1, 2, b = 1, 2) (28) 

holds. 

To construct a state with a fixed wave vector k, we first define translation 
operators T a in the directions u a (a = 1,2). For any state | $) the translation 
operators T\ and T 2 produce new states if | 4>) is not an eigenstate of them. We 
define, with n a = 1, 2, • • • (a = 1, 2), 

| $ (ni>n2) ) = T^TT I $> ■ (29) 
Using these states we make an eigenstate common to Ti and T 2 , 

| =CJ2\ ^(n u n 2 ))e- fk - {nA+nA) (30) 

"1,12 

Here the sum runs over all n\ and n 2 for which | ^( rai! n,))'s are linearly indepen- 
dent. The normalization factor C in (|30l) is therefore \j\f~N when all | $( ni .„ 2 )) 
are independent to each other. It is easy to see that 

f a | *(fc)) =| *(fe))e <E - <t » (a = 1,2), (31) 
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with the periodic boundary conditions in both directions. 

In order to describe our sets of basis states we consider a state | s) =\ 
sf , s\, ■ ■ ■ , s z N ) having a fixed z component s\ = +1/2 or —1/2 on each site i of 
the lattice. Then 2 N | s)'s form a set of basis states {| s)}. 

The set of basis states {| j)} is a complete set constructed from | s)'s for 
which the z-component of the total spin S z = J2i s i — 0- Each state {| j)} 
is defined by a linear combination of several | s)'s so that it has translational 
symmetries for the zero momentum and even under the n rotation, 2tt/3 rotation 
and the reflection. Total number of the basis states, namely Ny, in this set 
amounts - N C N/2 /l2N - 5.6 x 10 10 for the N = 48 lattice^ 

The set of basis states {| j(k))} should be constructed for each wave vector 
k. A basis state | j(k)) is determined by a state | s) in the S z = 1 sector, using 
it as | $) in (|2"5|) . Then the set of basis states {| j(k))}, which fulfil 

f a \j(k))=\j(k))e^ (a = 1,2), (32) 

is completed with as many states | s)'s as to span the full Hilbert space with 
S z = 1- Note that with this set of basis states one should use complex coefficients 
in expansions of states in the calculations. The total number of the basis states 

Z N CV/2-1/ 



for the 48-site lattice isH N v (k) = N C N/2 -i/N - 6.5 x 10 11 
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Figure 1: Magnetization on the sublattice defined by yM/v ■ The filled di- 
amond is our result on the 48-site lattice which we estimate by 50 samples of 



R 



(m) 
MO 



in (|T7|) with e = 4 x 1CT 3 and m — 15, using the operator O in (fH?)) . The 
statistical error is much smaller than the mark. Data on smaller lattices [ID] and 
results from the LSWT[4"T]. dashed line by (j21~j) . are also plotted. 




Figure 2: Brillouin zone on the triangular lattice and the wave vectors of the 
48-sitc lattice. 
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k/kO 

Figure 3: Static structure function, (S(—k) ■ S(k)) normalized by 8/(N + 6) for 
iV-site lattice. We obtain data on a 48-site lattice (filled diamonds) by R^o in 
p7|) with e = 5 x 10~ 2 and m — 15. Here the opertor O is given by (p?!?|) and 
fco = 47r/3. The number of samples is 50, with which statistical errors are quite 
small compared to marks. Squares are data for the 36-site system from ref.[6]. 
Solid lines are to guide eyes. 
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Figure 4: Magnon spectra along a path OQDBO in Fig. [21 Crosses are our data 
A(k X7 k y ) estimated on the 48-site lattice by one magnon gap energies ([25]) . 
Statistical errors and a solid line to guide eyes are also presented in the figure. 
The dashed line is the result from the LSWT[20]. 
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